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Abstract 



We examine the strong coupling limit of both compact and non-compact 
quantum electrodynamics (QED) on a lattice with staggered Fermions. We 
show that every SU(./Vl) quantum antiferromagnet with spins in a particu- 
lar fundamental representation of the SU(A ? l) Lie algebra and with nearest 
neighbor couplings on a bipartite lattice is exactly equivalent to the infi- 
nite coupling limit of lattice QED with the number of flavors of electrons 
related to N L and the dimension of spacetime, D + 1. There are N L 2- 
component Fermions in D = 1, 2A^l 2-component Fermions in D = 2 and 
27Vl 4-component Fermions in D = 3. We find that, for both compact and 
non-compact QED, when N L is odd the ground state of the strong coupling 
limit breaks chiral symmetry in any dimensions and for any iV L and the con- 
densate is an isoscalar mass operator. When N-^ is even, chiral symmetry is 
broken if D > 2 and if is small enough and the order parameter is an 
isovector mass operator. We also find the exact ground state of the lattice 
Coulomb gas as well as a variety of related lattice statistical systems with 
long-ranged interactions. 



1 Introduction 



It is reasonably well established that, as the bare coupling constant of mass- 
less quantum electrodynamics (QED) is increased, there is a phase transition 
which breaks chiral symmetry and generates an electron mass. The mech- 
anism is similar to that in the Nambu-Jona-Lasinio model IJ where chiral 
symmetry breaking occurs when the four-Fermion interaction is sufficiently 
strong and attractive that a bound state and the resulting symmetry break- 
ing condensate forms. In the case of QED it is the Coulomb attraction of 
electrons and positrons which, as the electric charge is increased, gets strong 
enough to form a condensate. In QED, this phase transition has been seen 
in both three and four spacetime dimensions using numerical as well as ap- 
proximate analytic techniques. In four dimensions, numerical simulations of 
compact and noncompact || lattice QED with dynamical Fermions indi- 
cate presence of a phase transition as the bare electric coupling is increased. 
For compact QED the transition is first order and for non-compact QED it 
appears to be of second order. 

For non-compact QED this phase transition is also found in the continuum 
using approximate analytical techniques such as the solution of Schwinger- 
Dyson equations in the quenched ladder approximation ||, ^ |£| |7], RL |?| . The 
critical behavior has the additional interesting feature that certain pertur- 
batively non-renormalizable operators such as four-Fermion operators can 
become relevant there [|[ [|. 

In three spacetime dimensions, a similar behavior is found in the large iV 
expansion ITD, ITTJ where the inverse of the number of Fermion flavors, 1/N, 



plays the role of coupling constant. Both analytical techniques [[T(J, |TTJ and 
numerical simulations [12] find a critical value of N, above which QED is 



chirally symmetric and below which the chiral symmetry is broken. 

Furthermore, there are some recent analytic proofs that the strong cou- 
pling limit of QED breaks chiral symmetry. Salmhofer and Seiler [13] showed 
using a Euclidean spacetime approach and staggered Fermions that in four 
dimensions, four-flavor QED (as well as some other U(iVc) gauge theories 
with iVc < 4) has a chiral symmetry breaking ground state when the electric 
charge is infinite. Subsequently, using staggered Fermions in the Hamiltonian 
approach it was shown [14, 15, 16] that the strong coupling limit of QED has 
a chiral symmetry breaking ground state in any spacetime dimension greater 
than two and when there are a specific number of Fermion flavors (two four- 



1 



component in 2+1 and four four- component Fermions in 3+1 dimensions, 
also for other U(iVc) gauge theories with any iVc). This is an elaboration of 
previous arguments for chiral symmetry breaking [TB| using strong coupling 



Hamiltonian methods (for a review, see |19[|). Also, it is partially motivated 



by an interesting previous work of Smit |2(| where he uses naive and Wilson 
Fermions (and also finds a particular kind of mapping to an SU(iV) antifer- 
romagent) to analyze chiral symmetry breaking in QED. 

If the chiral symmetry breaking phase transition in four dimensions is 
of second order, as it seems to be in non-compact QED ||, it provides 
a nonperturbative zero of the beta-function for the renormalized coupling 
constant. Furthermore, if the critical behavior there differs from mean field 
theory, the resulting ultraviolet fixed point would allow QED, in the limit 
of infinite cutoff, to avoid the Landau pole, or Moscow zero [[H], |22|, ^3 



which otherwise renders it trivial. It is an interesting and nontrivial question 
whether, via this mechanism, QED could be an example of a nontrivial field 
theory which exists in four dimensions. 

The existence of the strong-coupling phase transition has long been advo- 
cated by Miransky |5|, |§. He visualized the mechanism for chiral symmetry 
breaking as a "collapse" of the electron-positron wave-function, similar to the 
behavior of the supercritical hydrogen atom with bare proton charge Z>137. 
He argued that there are two ways to screen supercritical charges. In the case 
of a supercritical nucleus, the high electric field produces electron-positron 
pairs, ejects the positrons and absorbs the electrons to screen its charge. In 
the case of supercritical electron and positron charge, he argued that the 
pair production is suppressed by Fermion mass, so a system can stabilize 
itself by increasing the electron and positron masses; thus the tendency to 
break chiral symmetry. Miransky studied continuum QED in the quenched 
approximation, using the ladder Schwinger-Dyson equations to sum the pla- 
nar photon-exchange graphs. In this approximation, there is a line of critical 
points, beginning at bare coupling e 2 = and ending at e 2 c = tt/3 where the 
theory breaks chiral symmetry dynamically and has an interesting contin- 
uum limit. A behavior very similar to this was found in lattice simulations 



using quenched Fermions[£J, |25| . At least some of its qualitative features are 
expected to survive the presence of Fermion loops in realistic QED. 

In this Paper we shall examine Miransky' s collapse phenomenon in a 
more physical context by studying a lattice version of QED which is similar 
to a condensed matter system. We do this by arguing that the lattice ap- 
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proximation to QED with staggered Fermions in any number of dimensions, 
D > 1, resembles a condensed matter system of lattice electrons in an exter- 
nal magnetic field and with a half-filled band (this argument was also given 
in [14], [H| 0). In that picture, the breaking of chiral symmetry and the ana- 
log of the collapse phenomenon is the formation of either charge or isospin 
density waves and the resulting reduction of the lattice translation symmetry 
from translations by one site to translations by two sites. It occurs when the 
exchange interaction of the electrons, which is attractive, dominates the ten- 
dency of the direct Coulomb interaction and the kinetic energy to delocalize 
charge, giving an instability to the formation of commensurate charge den- 
sity waves. This forms a gap in the Fermion spectrum and a particular mass 
operator obtains a vacuum expectation value. This gives an intuitive picture 
of how strong attractive interactions in a field theory can form a coherent 
structure. Here, the commensurate density waves in the condensed matter 
system correspond to a modulation of the vacuum charge or isospin density 
at the ultraviolet cutoff wavelength in the field theory. 

We are interested in massless quantum electrodynamics with action 

S = j d D+1 x [-—F^ + E^VK + J (1) 

where A is the ultraviolet cutoff, e is the dimensionless electric charge and 
there are Np flavors of 2^' D+1 ) //2 l-dimensional Dirac spinors. (Here, [(Z)+l)/2] 
is the largest integer less than or equal to (D + l)/2.) In even dimensions 
(when D+l is even), the flavor symmetry is SUr(Nf) x SUl(Nf). In odd 
dimensions there is no chirality and the flavor symmetry is SU(Np). What 
is usually referred to as chiral symmetry there is a subgroup of the flavor 
group. We shall use a lattice regularization of flip and study the limit as e 2 — ► 
oo. We shall use Nl flavors of staggered fermions. In the naive continuum 
limit, this gives Np = Nl 2-component fermions in D=l, Np = 2Nl 
2-component fermions in D=2 and Np = 2Nl 4-component fermions in 
D=3. Though the lattice theory reduces to ([I]) in the naive continuum limit, 
the lattice regularization breaks part of the flavor symmetry. All of the 
operators which are not symmetric are irrelevant and vanish as the lattice 
spacing is taken to zero. In D=3, the SUr(Nf) x SUl(Nf) symmetry is 
reduced to SU(Np/2) (in D=3, Nl = Np/2) and translation by one site 
in the 3 lattice directions. These discrete transformations correspond to 



3 



discrete chiral transformations in the continuum theory. In D=2, the lattice 
has SU(N F /2) (in D=2, N L = N F /2) symmetry and two discrete ("chiral") 
symmetries. In D=l, the SU L (N F ) x SU R (N F ) symmetry is reduced to 
SU(N F ) (in D=l, N F = N F ) and a discrete chiral transformation. In each 
case, the discrete chiral symmetry is enough to forbid Fermion mass and it 
is the spontaneous breaking of this symmetry which we shall examine and 
which we shall call "chiral symmetry breaking". We shall also consider the 
possibility of spontaneous breaking of the SU(N F ) flavor symmetry. In the 
continuum the SU (N F ) corresponds to a vector-like symmetry subgroup of 
the full flavor group. Since the lattice and continuum theories have different 
symmetries, the spectrum of Goldstone bosons, etc. would be different in 
the two cases. In the following, we shall not address the source of these 
differences but will define QED by its lattice regularization and discuss the 
realization of the symmetries of that theory only. 

In this Paper we shall prove that, in the strong coupling limit, lattice 
QED with Al lattice flavors of staggered Fermions is exactly equivalent 
to an SU(Al) quantum antiferromagnet where the spins are in a particular 
fundemental representation of the SU(./Vl) Lie algebra. Furthermore, mass 
operators of QED correspond to staggered charge and isospin density oper- 
ators in the antiferromagnet. Thus, the formation of charge-density waves 
corresponds to chiral symmetry breaking and the dynamical generation of 
an iso-scalar Fermion mass whereas Neel order of the antiferromagnet corre- 
sponds to dynamical generation of Fermion mass with an iso-vector conden- 
sate. As a result of this correspondence we can use some of the properties of 
the quantum antiferromagnets to deduce features of strongly coupled QED. 

We find that, in the infinite coupling limit, the properties of the electronic 
ground state of compact and non-compact QED are identical. The fact 
that compact QED confines and non-compact QED (at least in high enough 
dimensions) does not confine affects only the properties of the gauge field 
wavefunctions and the elementary excitations. 

We find that the case of even Al and odd A^ are very different: 

When Al is odd, the vacuum energy in the strong coupling limit is pro- 
protional to e 2 , the square of the electromagnetic coupling constant. We 
also find that chiral symmetry is broken in the strong coupling limit for all 
odd A" L > 1 and for all spacetime dimensions D + 1 > 2. The mass operator 
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which obtains a nonzero expectation value is a Lorentz and flavor Lie algebra 
scalar. There are also mass operators which are Lorentz scalars and which 
transform nontrivially under the flavor group which could get an expectation 
value and break the flavor symmetry spontaneously if is small enough. 
We expect (but do not prove) that when N L increases to some critical value 
there is a phase transition to a disordered phase. 

In contrast, when is even we find that the vacuum energy in the 
strong coupling limit is of order 1. We find that chiral symmetry is broken 
when the spacetime dimension^ D + 1 > 3 and when Ni, is small enough. 
The mass operator which gets an expectation value is a Lorentz scalar and 
transforms nontrivially under the flavor SU(./Vl). Thus, flavor symmetry is 
spontaneously broken. As in the case of odd N-^, we expect that there is an 
upper critical where there is a transition to a disordered phase. 

In the course of our analysis, we find the exact ground state of the gen- 
eralized classical Coulomb gas model in D-dimensions with Hamiltonian 

e 2 

#coui = J! p( x )9(% - y)p(y) (2) 

<x,y> 

where the variable p(x) lives at the sites of a square lattice with spacing one 
(with coordinates (xx, . . . , Xd) and integers) and takes on values 

Nl Nl j_ 1 Nl ^1 Nl AT AA ■ + ^ 

- -y, — 2" + ' " " " ' ~2~ ~2~ an mte S er ( 3 ) 

and the interaction is long-ranged 

g(x — y) ~ \x — y\ 2 ~ D , as|x — y\ — > oo (4) 

In any space dimension, D > 1 we prove that there are two degenerate ground 
states which have the Wigner lattice configurations 

p(z) = ±I(-l)£f=^ (5) 

When iVx = 1 this is a long-ranged Ising model which is known to have 
antiferromagnetic order, even in one dimension. It is also in agreement with 
a known result about the Ising model in two dimensions [^(| . 

4 1+1 dimensions is a special case which we will discuss later (Section 4). 
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The difference between even and odd A^ can be seen to arise from a cer- 
tain frustration encountered when, with odd Ni, one simultaneously imposes 
the conditions on lattice QED which should lead to charge conjugation invari- 
ance, translation invariance and Lorentz invariance of the continuum limit. 
This frustration is absent when Al is even. This is not an anomaly in the 
conventional sense of the axial anomaly or a discrete anomaly encountered in 
the quantization of gauge theories, as no symmetries are incompatible, but 
it is nevertheless an interesting analog of the anomaly phenomenon. 

Note that in 3 spacetime dimensions the difference between even and 
odd Al is the difference between an even and odd number of 4-component 
Fermions. Our result that for an odd number of 4-component Fermions the 
chiral symmetry is always broken for large coupling seems to contradict the 
continuum analysis in [pi], [TT[1 . We do not fully understand the reason for 
this, but speculate that it is related to the lattice regularization. Their model 
is very similar to ours in that they effectively work in the strong coupling 
limit when they replace the ultraviolet regularization which comes from hav- 
ing a Maxwell term in the QED action by a large momentum cutoff. In our 
case we have a lattice cut-off and the result should be very similar. Note 
that we agree with them when Al is even and there is an even number of 4- 
component Fermions. In that case, we are also in qualitative agreement with 



recent numerical work |I2J which, since it uses Euclidean staggered Fermions 
and there is a further Fermion doubling due to discretization of time, can only 
consider the case where there is an even number of 4-component Fermions. 
The anomalous behavior that we find with odd numbers of 4-component 
Fermions exposes a difficulty with treating the flavor number iV as a contin- 
uous parameter. 

A hint as to why 2+1-dimensional Fermions should come in four-component 
units appears when we formulate compact QED in the continuum using the 
SO (3) Georgi-Glashow model with spontaneous breaking of the global sym- 
metry, 50(3) — > U(l) and the limit of large Higgs mass [27]. We begin with 
the model 

£ = A E(^) 2 + \ E wr ■ wr + t((E w - v2 f ( 6 ) 

4e a=l / a=l 4 a=l 

The spectrum contains a massless photon which in the A — > oo and v 2 — > oo 
is the only light excitation. Since the U(l) gauge group is a subgroup of 
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S0(3), it is compact. We wish to make 2+1-dimensional electrodynamics by 
coupling this model to Fermions. It is known that, if the resulting theory 
is to preserve parity and gauge invariance simultaneously, we must use an 
even number of two-component Fermions which, in the minimal case are also 
SU(2) doublets [fH 



£ F = i/}(i>y- V + 7- A + ^(m + g<p- a))^ (7) 

Here, the Fermions are four-component and have a parity invariant mass and 
Higgs coupling. Each SU(2) doublet contains two electrons (which can be 
defined so that both components of the doublet have the same sign of electric 
charge), thus the basic Fermion for compact QED has eight components. The 
maximal chiral symmetry is obtained in the massless case with vanishing 
Higgs coupling. It is possible, by suitable choice of mass and Higgs coupling 
i.e. m = ±g < <fi >, (and also reduction of the chiral symmetry) to make four 
of the Fermions heavy, leaving four massless components. This is consistent 
with parity and gauge invariance. It is also interesting to note that staggered 
Fermions on a Euclidean lattice, where time is also discretized, produce eight 
component continuum Fermions. 

In the lattice gauge theory, the coupling of gauge fields to the electron field 
does not distinguish between compact and non- compact QED. Therefore, 
the above considerations for the Fermion content should apply to both cases. 

Similarly, in four dimensions it is known that the Fermion multiplicity ob- 
tained by lattice regularization is the one which is compatible with the axial 
anomaly p9[ |3(J. An alternative way to see this constraint is to again form 
compact QED using the Georgi-Glashow model. Then, because of Witten's 
SU(2) anomaly JHj], we are required to use at least two (and generally an even 
number of) 2-component Weyl Fermions which are at least SU(2) doublets. 
This seems to indicate that the minimal Fermion would have eight compo- 
nents. This is the minimal number that we find using staggered Fermions 
in the Hamiltonian approach (see Appendix B). It is interesting that in Eu- 
clidean staggered Fermions produce sixteen components, and therefore in our 
terminology always lead to the case with even N^. 

A different context where a correspondence between spin systems and 
gauge theory appears is in the study of strongly correlated electron systems 
in condensed matter physics. An important issue there is the correspondence 
of spin systems such as the quantum Heisenberg antiferromagnet (which is 
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equivalent to the strong coupling limit of the Hubbard model at 1/2 filling) 
and certain limits of various kinds of lattice gauge theories. Gauge theory- 
like states can be obtained as an approximate low energy theory in the mean 
field approximation. An example is the "flux phases" suggested in |32], |33], ^4 



and also various dimerized phases J35] which are disordered and under some 
conditions compete with the Neel phase, particularly in 2+1 dimensions (for 
a review, see |[(J). Affleck and Marston J34J showed how to get the flux phase 
from mean field theory. It is a locally stable but not global minimum of the 
free energy of an antiferromagnetic spin system and could presumably be 
stabilized by adding certain operators to the Hamiltonian. The low-energy 
limit in this phase resembles strongly coupled 2+1-dimensional lattice QED 
with four species of massless 2-component electrons (because of the two spin 
states of the lattice electrons, this is the case N L = 2 in 2+1-dimensions). 

The picture that we shall advocate in this paper is that there should be 
a phase transition between the flux phase and the Neel ordered phase of the 
antiferromagnet. This phase transition is governed by the strength of the ef- 
fective electromagnetic coupling constant. The Heisenberg antiferromagnet, 
which is known to have a Neel ordered ground state, is obtained in the limit 
of infinite electric coupling constant. For weaker couplings, the system can 
be in the flux phase where the electrons are massless. As the coupling is 
increased, the Neel state is recovered by formation of a commensurate spin 
density wavef], which corresponds to spontaneous chiral symmetry breaking 
and the generation of iso-vector Fermion mass in the effective QED. (Spin 
in the condensed matter system corresponds to iso-spin in effective QED.) 

In Section II we discuss the formulation of QED on a lattice. We give 
a detailed discussion of discrete symmetries and also of gauge fixing which 
is necessary to make non-compact QED well-defined. In Section III we 
discuss the strong coupling limit and show how quantum antiferromagnets 
are obtained in the strong coupling limit for both compact and non-compact 
QED. We discuss the properties of the electron ground state for both cases of 
N L even and odd. We also discuss the implications of the mapping between 
antiferromagnets and strong coupling gauge theories for the symmetries of 
the ground state. Section IV is devoted to concluding remarks. 

We review some of our notation in Appendix A and the essential features 

5 This is a commensurate charge density wave for each spin state whose phase is such 
that the condensate has zero electric charge. 
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of the staggered Fermion formalism, with emphasis on those aspects which 
are important for our arguments, in Appendix B. Appendix C is devoted to 
a review of the Fermion formulation of spin systems. 
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2 QED on a Lattice 



In this Section we shall set the Hamiltonian formalism of Abelian gauge fields 
on a lattice. For the most part, this formalism can be found in some of the 
classic reviews of lattice gauge theory, |L9| for example. A novel feature 
of the present Section is a careful treatment of the normal ordering of the 
charge operator and a discussion of the ensuing discrete symmetries. This 
normal ordering turns out to be important if the continuum limit is to have 
the correct behavior under C, P and T transformations. It will also be 
important in our later solution of the strong coupling limit. 

2.1 Hamiltonian and Gauge Constraint 

We shall discretize space as a cubic lattice and, in order to use the Hamil- 
tonian formalism, time is left continuous. We use units in which the lattice 
spacing, the speed of light and Planck's constant are all equal to one. (See 
Appendix A for a summary of our notation and Appendix B for a review 
of staggered Fermions.) Lattice gauge fields are introduced through the link 
operators 

Ui(x) = e iA ^ x) (8) 

which correspond to the link [x, i] of the lattice. Electric fields propagate on 
links of the lattice and the electric field operator Ei(x) associated with the 
link [x,i] is the canonical conjugate of the gauge field 

[A i {x) i E j {y)] = iS ij 6{x-y) (9) 

The gauge field and electric field operators obey the relations 

A^(x) = -Ai(x - i) , E^{x) = -Ei(x - i) (10) 

The Hamiltonian (of non-compact QED) is 

[x,i] [x,i,j\ 

+ E (W 1 ^ + i)e iA ^r{x) + h-c) (11) 

[x,i] 
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where the second term contains a sum over plaquettes and the magnetic field 
is defined as the curvature of the gauge field at the plaquette [x,i,j], 

B[x,i,j] = Ai{x) + Aj{x + %)+ A-i{x + A-j{x + 3) 

= Ai{x) - Aj{x) + Aj(x + i) - Ai(x + j) (12) 

As is discussed in Appendix B, the hopping parameter t[ x $ contains phases 
which produce a background magnetic flux it per plaquette. In the weak 
coupling continuum limit, the magnitude of t\ X £ is one, \t\ 2 = \t\ x {\\ 2 = 
1 in order that the speed of the free photon and free electron fields are 
equal, i.e. so that the low frequency dispersion relations for both the photon 
and electron have the same speed of light. However, in order to obtain a 
relativistic continuum limit in general it is necessary to make \t\ a function 
of e 2 . We shall find that in the limit where e 2 is large, the speed of light is 
proportional to \t\/e and it is necessary that \t\ ~ e as e 2 — > 00. 

The Hamiltonian we have written in (|TT|) is appropriate to non-compact 
QED. If we wish to study compact QED we must make the Hamiltonian 
symmetric under the field translation 

Mx) -> Ai(x) + 2tt (13) 

for any [x,i]. This is accomplished by replacing the second term in flTTD by 
a periodic function of the magnetic flux, so that (for compact QED) 

H c = J2jE t (xy+ £ |sin 2 (£[*,z,j]/2) 

[x,i] [x,i,j] 

+ E ( W f (z + i)e iA ^H a (x) + h.c.) (14) 

[x,i] 

Both (|ll|) and (0) reduce to the standard Hamiltonian of QED in the naive 
weak coupling continuum limit. Away from that limit the behavior of the 
dynamical systems described by the two Hamiltonians can differ significantly. 
For example, in the strong coupling limit compact QED is a confining p7| 



theory whereas non-compact QED in not confining. Also, the phase transi- 
tion seen in numerical simulation of the two theories differ. In the compact 
case the phase transition associated with chiral symmetry breaking is gen- 
erally of first order whereas it is second order for the non-compact case 
[HI. The source of some of these differences generally have to do with the 
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symmetry (|TtJ). In ( |TT| ) and ( PD we have introduced 7V"l flavors of lattice 
Fermions labelled by the index a = 1, . . . , N-^. 

In both compact and non-compact QED, the Hamiltonian is supple- 
mented with the constraint of gauge invariance. The gauge transformations 
of the dynamical variables, 



Ai(x) + ViA(x) 
Ei(x) 



A : 


A(x) 


A : 


Ei(x) 


A : 




A : 





^(x)e- iA{x) 
are generated by the operator 

Q A = Y1 (-ViA(x)Ei(x) + A(x) (r\x)r(x) - N L /2 



(15) 



(16) 



The local generator of gauge transformations where A has compact support 
is 



eg. 



A 



dA{x) 



G(x) = V • E{x) + ^{x'W a {x) - N L /2 



(17) 



Both ([TBI) and (|i~7[ ) commute with the Hamiltonians in ( |TT1) and ([14|) . 

In ([TB]) and flT7|) we have subtracted the constant from the charge 

density operator in order to make the gauge generator odd under the usual 
charge conjugation transformation 



A(x) 
Ei(x) 

tp a (x) 



- A(x) 

- Ei(x) 

(-i)Ef=i^^( x ) 



(18) 



In fact, the Fermionic charge term in ( JTBj ) can be put in the manifestly odd 
form ^[i[) a <(x), ip a {x)} Of course, charge conjugation symmetry of the lattice 
theory is necessary to obtain charge conjugation of the continuum theory. 
We shall see later that, particularly at strong coupling, the subtraction term 
in ([TBI) plays an important role. It seems to have been ignored in previous 
literature ( for example, see ]18|, |19[). Its presence is particularly important 
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when Ni, is odd since the charge operator has no zero eigenvalues in that 
case (the eigenvalues of ift'if) are integers). 

Chiral symmetry is related to translation invariance by one site. The 



11) and 


(HD are 


invariant under the transformations 


Xj ■ 


A(x) - 


-> Ai{x + j) 


Xj ■ 


Ei(x) - 


-> Ei(x + j) 


Xj ■ 


^ a {x) - 


-> (-l^w^O^ + j) 


Xj : 




-> (-l)^^+ lXfc ^ af (x + j) 



(19) 

for j = 1, • • • , D. 

In the following we shall use the charge conjugation symmetry which is a 
combination of these two transformations: 

C = txi (20) 

This is necessary if the mass operators which we define in Appendix B is to 
be invariant under charge conjugation symmetry. Also, we shall see that the 
strong coupling ground state is invariant under C but not under either £ or 
Xj by themselves. 

The dynamical problem of Hamiltonian lattice gauge theory is to find 
the eigenstates of the Hamiltonian operator (|TT| ) or (|14]) and out of those 
eigenstates to find the ones which are gauge invariant, i.e. which obey the 
physical state condition (or, the "Gauss' law" constraint) 

£(a;)|^ phys . >=0 (21) 

Note that the gauge constraint and physical state condition are the same for 
both compact and non-compact QED. In the case of compact QED there 
is the additional symmetry (|T3| ) which, being a large gauge symmetry, can 
be represented projectively. In fact, when we later work in the Schrodinger 
picture we shall require that the quantum states which are functions of a 
configuration of the gauge field transform as 

\Ai(x) + 2ixn >= exp (in9[x, i)) \Ai(x) > (22) 

There is a separate parameter 9[x, i] for each link of the lattice. These pa- 
rameters originate in a way similar to the theta-angle in ordinary QCD. The 
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symmetry ( ]22j) together with the commutator (0) imply that the spectrum 
of the electric field operator has eigenvalues which are separated by integers 
and offset by 9: 

spectrum [E'j(x)] = {9[x, i] + integers} (23) 

The Hamiltonian and gauge constraints can be obtained from the gauge 
invariant Lagrangian 

L = Y,r ] {x){id t - A (x))r(x) + E Ei(x)Mx) 
x l x A 

+ £ Ei(x)ViMx) + E Mx)Nl/2 - H (24) 

[x,i] x 

where the temporal component of the gauge field has been introduced to 
enforce the gauge constraint and the time derivative terms give the correct 
symplectic structure. Note that, in order to get Lorentz invariance of the 
Fermion spectrum in the weak coupling (naive) continuum limit, we require 
half-filling of the Fermionic states, i.e. that the total charge defined by 



has zero vacuum expectation value. 

Here we have considered massless QED. As well as the gauge invariance 
and charge conjugation invariance discussed above, the Hamiltonian is sym- 
metric under the discrete chiral transformations (19) which on the lattice 



corresponds to symmetry under translation by one site. In later Sections, we 
shall consider the possibility of spontaneous breaking of this symmetry. 

2.2 Gauge fixing and quantization 

We shall quantize the gauge fields in the Schodinger picture. The quan- 
tum states are functions of the link operators which are represented by the 
classical variables Ai(x) and the electric field operators are derivatives 

£ - w s -ah) (25) 
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We shall also consider the usual Fock representation of the Fermion anticom- 
mutator. The empty vacuum is the cyclic vector 

^ a (x)\0>=0 Va,x (26) 

and Fermions occupying lattice sites are created by ip a ^(x). 

2.2.1 Compact QED 

In compact QED the spectrum of the gauge generator is discrete and a state 
which obeys the physical state condition can be normalized, thus imply- 
ing that there is no need for additional gauge conditions. The basis wave- 
functions for compact QED (in the basis where the Fermions density and 
eletric field operators are diagonal) are \l/[n(x)]$[y4] with the Fermion states 

V[n(x)] = II II (x)) na(x) |0 > (27) 

X 0=1 

labeled by vectors n(x) = (ni(x), • • - , 71jv l (:e)) with n a (x) = or 1, and the 
photon states 

V[A}=exp(tJ2e t (x)A t (x)) (28) 

where the eigenvalues ei(x) of the electric field operator are in spectrumf-E^x)] 
(|23|). Furthermore the states of the photon field are normalized using the in- 
ner product 

< d>xLA],<f> 2 L4] >= I] I*" ^&9\[A]* 2 [A] (29) 
; ., jo 1~ 

[x,i] 

and the Fermion states have conventional inner product given by < 1 >= 1 
and the canonical anticommutator relations of the fermion field operators. 
The physical state condition fl2T|) gives the additional restriction that 

V^x) = -p(x) = -E EKW - 1/2) (30) 

x a 

where p(x) is the charge density (i.e. the eigenvalue of ip a ' (x)ip a (x) — N-^/2). 
Pictorially, we can think of this as containing lines of electric flux joining sites 
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whose charges are non-zero and also closed loops of electric flux. In the strong 
coupling limit, the Hamiltonian is diagonal in the basis (p8|) . This gives a 
natural explanation of confinement in compact QED in the strong coupling 
region. If we add a particle-antiparticle pair to a state in ( p8|) it must be 
accompanied by at least a single line of electric flux. The energy of such 
a line of flux is proportional to its length. Therefore the electron-positron 
interaction grows linearly with distance and is confining. This is in contrast 
to the situation in non-compact QED where the electric flux is not quantized. 
In that particle-antiparticle pair can have many lines with arbitrarily 

small flux. The energy of the field is minimized by the usual Coulomb dipole 
configuration. In high enough dimensions this is not a confining interaction. 



2.2.2 Non-Compact QED 

In contrast to the case of compact QED, in non-compact QED the generator 
Q{x) fll7D of gauge transformations has a continuous spectrum. In order 
to obtain a normalizable ground state it is therefore necessary to fix the 
additional gauge freedom. 

In order to separate the gauge orbits of the photon field, we shall need to 
define the transverse and longitudinal components of the gauge fields. The 
transverse projection operator is 

Tij = <5y - ^ (31) 
V ■ V 

and the longitudinal projection operator is 

U 3 = (32) 
V-V 

They have the usual property of projection operators, 

T 2 = T L 2 = L TL = = LT (33) 

and also, 

1 = T + L (34) 

The transverse and longitudinal parts of the electric and gauge fields are 
obtained by 

Aj(x) = TijAj(x) Af(x) = LijAj(x) (35) 
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Ej{x) = TijEjix) Ef(x) ee /.,,/•-• t.r) (36) 

Note that there is an ambiguity in the precise definition of the projection 
operators T and L on a finite lattice due to the zero mode of the lattice 
laplacian. We have fixed this so that constant fields, Ai(x) = Oj = const., 
are purely transverse, <Zj = aj (see Appendix A). 

In order to quantize, it is necessary to solve the gauge constraint 

Q(x) = V-E{x)+p{x) ~0 (37) 

with 

p{x) ee ^ at (x)^ a (x) - iV L /2 (38) 



This is most easily accomplished by gauge fixing. The procedure |Jy[] is to 
find a gauge fixing condition which has a nonzero commutator with the gauge 
constraint. An example is the Coulomb gauge condition 

X{x) = V ■ A(x) ~ (39) 

The commutator of the gauge condition with the gauge generator is 

[ x ( x ),g( y )] = -iV.V5{x-y) (40) 

which is a non-degenerate matrix. We can then solve the two constraints by 
eliminating the longitudinal parts of the gauge and electric fields, 

A\{x) = (41) 



El{x) = -Vi— — -p(x) ee -^ Vi (x|— — f-\y)p(y) (42) 

V * V y V " V 

The remaining degrees of freedom obey canonical commutation relations 
which are derived from Dirac brackets. These brackets project the canonical 
Poisson brackets onto the constrained phase space. Given a set of constraints, 
£a with a non-degenerate Poisson bracket, det £b}pb 0; Dirac brackets 
for dynamical variables are obtained from Poisson brackets as 

{P, Q}db = {P, Q}pb ~ {P, U} PB iU Zb}pb {Zb, Q} pb (43) 
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In the present case, the brackets of the remaining variables are not modified. 
The commutator for the transverse photon and electric fields is 

[Aj{x),E]{y)]=iT ij 8{x-y) (44) 

The Hamiltonian depends only on the transverse photon and electric field 
and on the charged Fermion fields, 

# = E e ^{Ejix)f + E 2^* T (*)(-v • v)A T (*) 

[x,i] [x,i] 

+ E *M + i)e U7(x) r(x) + h.c.) (45) 

[x,i] 

The Coulomb interaction has appeared as a result of the solution of Gauss' 
law. 



This gauge fixing, which we have done following Dirac's procedure [p7 



can always be implemented by a canonical transformation. Gauss' law is 
solved by taking the ansatz for the physical states 

|phys >= exp ( iJ2 M x ) ■ Vi l -^-p{x) J \A T >= U\A T > (46) 

The exponential operator in fl46|) generates the unitary transformation 

Ai(x) = UAi{x)U ] = Ai(x) 
Ei(x) = UE^x)U ] = Ei(x) - V>J^p(x) 

tp(x) = Ui)(x)U ] = exp ^ 2 V ■ A(x)j ip(x) 

${x) = Uifi(x)lfl = ^{x) exp (-^—V ■ A{x)j (47) 

Note that the transformation of ip and ip* removes the longitudinal part of 
the gauge field from the covariant hopping term, 

${x + i)e iA * (x) ^{x) = ijj\x + i)e iA ^ x) tfj{x) (48) 
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Substituting the canonically transformed fields into the Hamiltonian yields 
the Hamiltonian fl45|) which is decoupled from A\{x) together with the canon- 
cially transformed Gauss' law which now states that E L = 0. The longitu- 
dinal parts E L and A L can now be dropped from the phase space and the 
resulting quantum theory is the same as that obtained from Dirac's procedure 
for solving the constraints. 

For D = 1 (1+1 dimensional spacetime), fl45f ) is the Hamiltonian of the 
lattice Schwinger model. It is worth pointing out that in that case, A T and 
E T just comprise one quantum mechanical degree of freedom, 

A T = J2A(x), E T = J2E(x) 

X X 

This corresponds to the fact that the only physical degree of freedom of the 
photon field on a 1 dimensional compact space is the Wilson loop (U(l) 
holonomy) exp(z J2 X A(x)). 

3 Strong Coupling Expansion 

Although the results are very similar, the implementation of the strong cou- 
pling expansion is somewhat different in the two cases of compact and non- 
compact QED. We shall treat the two cases separately. 

3.1 Non-Compact QED 

The conventional strong coupling, e 2 — > oo limit is difficult to implement for 
noncompact QED since the leading terms in e 2 in the Hamiltonian have a 
continuum spectrum. The alternative, but related procedure is the hopping 
parameter expansion, i.e. an expansion in the parameter \t\ in equation 
The terms in this expansion are very similar to a strong coupling expansion 
in that they contain inverse powers of e 2 . 

We begin by separating the Hamiltonian into two parts, a leading order 
part 

#o = E ^(Ej{x)f + £ ^ 2 Aj{x){-V ■ V)Aj(x) + 

[x,i] 

e 2 1 

EipW^I - ^ \y)p(y) ( 49 ) 

XM — V ' V 
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whose ground state we shall attempt to find exactly and a next-to-leading 
order part 

#i = E *M (^0* + i)e*WV°(aO + h.c.) (50) 
[«,»] 

which we treat perturbation. 

We first examine the structure of the ground state of H . First of all, it 
is a direct sum of the free transverse photon Hamiltonian and the Coulomb 
interaction which depends only on the Fermion operators. The wave-function 
therefore factorizes into a wavefunction for the free photon ground state and 
a wavefunction for the ground state of the four- Fermion operator in (ffi^). In 
the Schrodinger picture, the photon ground state is the Gaussian 



$photon[^] = 




(C the normalization constant) and the photon contribution to the ground 
state energy is just the ground state energy of the free photon theory and is 
of order zero in e 2 as well as \t\. 

The nature of the ground state of the four-Fermion part depends on the 
number of Fermion flavors, N^. In particular it is quite different when N-^ is 
even or odd and we shall treat these two cases separately. 



3.1.1 N L Even 

If A^l is an even number, the ground state of the operator 

e 2 1 

#coui = E -wp( x ) ( x \ —=—±:\y)p(y) ( 52 ) 

x,y Z — V ■ V 

is the state \g.s. > where 

N L 

E \g*. >= N L /2 \g.s. > (53) 

o=l 

i.e. with every site of the lattice half-occupied. It is easy to see that this is 
the case by noting that H cou i is a non-negative operator and that the states 
with zero charge density are zero eigenvalues of H con \. 
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This ground state is degenerate. At each site the quantum state is given 

by 



At/2 

J] ^ ait |0 > (54) 

i=l 

Since this quantity is antisymmetric in the indices oi, . . . , ajv L /2 it takes on 
any orientation of the vector in the representation of the flavor symmetry 
group SU(A^l) with Young Tableau made of one column with N L /2 boxes 
(m = N L /2 in Fig. 1). 



Figure l:The representation of SU(Nl) at each site when Nl is even. 

As in ref. |14], [L^, [L6| we observe that the degeneracy must be resolved by 
diagonalizing perturbations in the hopping parameter expansion. The first 
order perturbations vanish. The first non-trivial order is second order, 

6 2 = - < g.s.\H x H x \g.s. > (55) 
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This matrix element can be evaluated by noting that Hi creates an eigenstate 
of H different from the ground states with additional energy 



A£ = y + j(D)V 1 (x\-^\x) = e 2 (56) 

Diagonalizing the matrix of second order perturbations is equivalent to find- 
ing the spectrum of the effective Hamiltonian 



2\t 



2 



2\t 



[x,i] 



1 2 



e 2 

\X,l 



J2 J ab(x)Jba(x + i) (57) 



where the operators J a b(x) = ip a ^ (x)ip b (x) — |^ a fe, are the generators of the 
U(A^l) given in equation ( |C3D of Appendix C and obeying the Lie algebra in 
equation (|U2[). 



The constraint (p3|) on the total occupation number of each site, 

N L 
a=l 

reduces to SU(./Vl) (see Appendix C) and projects onto the irreducible repre- 
sentation given by the Young Tableau in Fig. 1. (This is one of the fundamen- 
tal representations of SU(ATl).) Furthermore, Q57D is just the Hamiltonian of 
the SU(./Vl) antiferromagnet in that representation. 

It is straightforward to see that the higher orders in the hopping param- 
eter expansion also have higher orders of 1/e 2 . In fact, if we consider the 
following limit, 

e 2 — > oo , \t\ 2 — > oo 

\t\ 2 /e 2 = constant (58) 

all higher order perturbative contributions to both the wavefunction and the 
energy vanish. Thus, in this limit, QED is exactly equivalent to an SU(A^l) 
antiferromagnet. That ( |58"D is the correct limit to take can be seen from the 
fact that, if the antiferromagnet in (|57|) is in an ordered state, the speed of 
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the spin-waves, which are the gapless low-energy excitations is proportional 
to \t\/e. They have linear dispersion relation uj{k) ~ \k\ and play the role of 
massless goldstone bosons for broken flavor symmetry. Their speed should be 
equal to the speed of light, which is one in our units. This implies that \t\/e 
should be adjusted so that the spin-wave spectrum is relativistic, u(k) = \k\. 
Hence the limit in (|j8|). 

When A/l = 2, this model is the quantum Heisenberg antiferromagnet in 
the j = 1/2 representation. It is known to have a Neel ordered ground state 
in D > 3 |H and there is good numerical evidence that it has Neel order in 



D=2. The antiferromagnetic order parameter is the mass operator 

E = E(" 1 )^= 1 Xi ^{x)ai/)(x) (59) 

X 

which obtains a vacuum expectation value in the infinite volume limit. Thus, 
when Nl = 2 the strong coupling limit breaks chiral symmetry and generates 
electron mass. It is interesting that in this case there is an iso-vector con- 
densate. In the strong coupling limit this seems unavoidable. The only way 
to get an iso-scalar condensate is with a charge density wave. However such 
a configuration always has infinite coulomb energy compared to an electric 
charge neutral but isospin carrying condensate. 

The low energy excitations of this systems (with energies of order \t\ 2 /e 2 
are spin waves. All other excitations have energies which go to infinity in 
the limit (|58|). The spin waves are the pions which are the scalar Gold- 
stone bosons arising from spontaneous breaking of the vector flavor symmetry 
SU{2) -+ U(l). 

For large there is some evidence that the SU(A r i J ) antiferromagnet in 
these specific representations has a disordered ground state [M. Particu- 
larly in 2+1-dimensions it is known that for infinite Nl the ground state is 



disordered |34[ Although it is beyond the scope of this paper, it would be 
intersting to investigate the A^ dependence of the ground state further. We 
shall comment on this in Section IV. 



3.1.2 N L Odd 

When N L is odd the charge density operator p(x) fl3lp which enters the 
Coulomb Hamiltonian (|52"D has no zero eigenvalues. Therefore the Coulomb 
energy of the ground state is necessarily of order e 2 for large e 2 . 
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Since the Coulomb Hamiltonian commutes with the charge density op- 
erator, p(x), they can be diagonalized simultaneously. The ground state of 
the Coulomb Hamiltonian should therefore also be an eigenstate of p(x). 
Therefore, to find the spectrum of (0) we consider all states which are also 
eigenfunctions of the local density, i.e. where at a given site x, 

p(x) = -N L /2, -N L /2 + 1, . . . , N L /2 (60) 

with the constraint of global neutrality 

E = 

X 

For convenience, we consider the system on a finite spatial lattice V = Vr 
with periodic boundary conditions (see Appendix A). Then the momenta k 
appearing in the Fourier transform are discrete. We consider the Coulomb 
Hamiltonian in momentum space 

e 2 1 

Z x,y — V ■ V 

e 2 1 1 

2 \V\ V 4££=ism (fei/2) 

where \V\ is the total number of spatial lattice sites and p(k) is the fourier 
transform of the charge operator p(x). Since < sin 2 (/cj/2) < 1 we can 
derive the lower bound for the Coulomb energy as 

Hcoul > —\p(k)\ 2 = —Yp(x) 2 > 

(the latter estimate follows from p(x) 2 > 1/4) This bound is saturated by 
the charge distributions 

\V\ 
~2 

where 7f = (tc, . . . ,ir) is the vector for which YliLi sin 2 (fcj/2) in the denomi- 
nator of (|BTD takes its maximum value. These are allowed configurations of 
the charge density, 

p (x) = ±J(-l)Sf=r 



Po{k) = ±— <%, 
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which give the ground state configurations of the Coulomb system. The 
electric field in these ground states is easily deduced from Gauss' law, 



These configurations break chiral symmetry in that they are not invariant 
under the transformation xi i n (0) but they are symmetric under charge 
conjugation C defined in Q2D|). 

The ground state energy per lattice site is 

£o(coui)/M = ^ (62) 

Note that it is of order e 2 . This is in contrast to the ground state energy when 
the number of lattice Fermion flavors is even, which is of order \t\ 2 /e 2 ~ 1. 

The ground states that we have found are highly degenerate in that only 
the number of Fermions at each site is fixed. Their quantum state can still 
take up any orientation in the vector space which carries the representation 
of the flavor SU(A^l) given by the Young Tableaux in Fig. 2. 
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Figure 2: Representation of SU(Nl) on each site of the even sublattice A 
and the odd sublattice B when Nl is odd. 

We have divided the lattice into two sublattices: A is all points where 
J2i %i is even and B where J2i x i is odd. Then, the differing occupation 
numbers on sites on each sublattice yield different representations of SU(Nl)- 

Again, this degeneracy must be resolved by diagonalizing the perturba- 
tions, which are non-zero in second order and the problem is equivalent to 
diagonalizing the antiferromagnet Hamiltonian (f)7]), this time with the rep- 
resentations depicted in Fig. 2. Also, in the limit (|58|) this correspondence is 
exact. 

The strong coupling ground states that we find when JVl is odd contains a 
charge density wave. The staggered charge density operator has expectation 
values 

< £(-1)^ V(aOV>0»0 >= ±1/2 (63) 

This condensate is an isoscalar and we have shown that it must always occur 
in all dimensions. When > 1 mass operators with certain generators of 
SU(A r L) could have expectation values if the ground state has antiferromag- 
netic order. However, unlike the case of even JVl, the antiferromagnetic order 
is not required in order to have chiral symmetry breaking. 

The ground state we find breaks chiral symmetry. This is a true dynami- 
cal symmetry breaking since, in infinite volume, the ground states which are 
related by a chiral transformation are never mixed in any order of strong cou- 
pling perturbation theory. Furthermore, there are no local operators which 
couple them. 

We conclude that the strong coupling ground state breaks chiral symme- 
try for any odd and in any dimensions. As in the case of even Nl there 
is also the possibility (and for small Nl the likelyhood) that the SU(N L ) 
antiferromagnet we obtain here is in a Neel state and the flavor symmetry is 
also broken. We shall not pursue this possibility here but refer the reader to 
the literature |3T| . 
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3.2 Compact QED 

The difference between compact and non-compact QED resides in the quan- 
tization of the gauge fields. In all cases the Fermionic state is identical in the 
two cases. In compact QED the eigenstates of the electric field operator are 
normalizable and can be used for the ground state. In this case we separate 
the Hamiltonian into three terms, 

[x,i] 

[x,i] 

H 2= E ^sm 2 (B[x,i,j)/2) (64) 

In the strong coupling limit it is necessary to solve Gauss' law for the 
electric fields and the charge distribution in such a way as to minimize H . 

When A/l is even, the charge operator has zero eigenvalues and Gauss' 
law has the solution where Ei(x) = 0, which is an obvious minimum of the 
Hq, and there are fermions on each site. This is similar to the situation 
in non-compact QED when ]Vl is even. Also, the degeneracy of this state 
must be resolved in the same way, resulting in the effective Hamiltonian 



(p7 ) which describes the SU(A/l) antiferromagnet in the representation with 
Young Tableau having one column with boxes shown in fig. 1. Again, 

we expect that this system has Neel order in D > 2 if A/l is small enough 
and the chiral symmetry of electrodynamics is broken, along with the SU(ATl) 
flavor symmetry. 

When A^l is odd, since the charge density operator has no non-zero eigen- 
values, it is impossible to find a zero eigenstate of the Gauss' law constraint 
operator without some electric field. The problem which we must solve is to 
minimize the energy functional E 2 subject to the constraint V • E — — p 
where at each site p has one of the values 

— — , — — + 1, . . . , — N L an odd integer 

It is straightforward to show that the charge distribution and electric field 
which one obtains is identical to those in the case of non-compact QED with 
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odd N L} 

Po(x) = i(-l)^- 4 E t {x) = -L(-1)EL** 

The ground state degeneracy is again resolved by diagonalizing perturbations 
and, again the true ground state of the strong coupling limit is the ground 
state of the effective Hamiltonian fl5T|) when the SU(./Vl) spins take on the 
configurations in Fig. 2. 

Notice that in the ground state, the electric fields are not integers, but 
on each link, the spectrum of the electric field operator is 1/4D+ integers. 
The "theta angles" 1/4D survive all orders in strong coupling perturbation 
thery. 



4 Remarks 

In this paper we have analyzed the possibility of chiral symmetry breaking in 
the strong coupling limit of quantum electrodynamics using the Hamiltonian 
picture and a lattice regularization. We chose to use staggered Fermions 
because they give the closest analog to interesting condensed matter physics 
systems. Also, unlike Wilson Fermions which, in the Hamiltonian picture, 
have no chiral symmetry at all, they have a discrete chiral invariance which 
forbids Fermion mass and it is sensible to ask questions about dynamical 
mass generation. 

In 1+1 dimensions, staggered Fermions give N-^ species of 2- component 
Dirac Fermions. When = 1 we obtain the Schwinger model with a lattice 
regularization. Also, in this case, we have found that the chiral symmetry is 
broken dynamically. Of course, due to the staggered Fermion regularization 
there is no continuous chiral symmetry, which is as it should be since it should 
be impossible to regularize the Schwinger model so that there is simultaneous 
continuous chiral and gauge symmetry. However, to match the solution of the 
continuum Schwinger model, the Fermion should obtain mass. This indeed 
happens in our strong coupling limit by spontaneous symmetry breaking. 
(Although we disagree with some aspects of the formalism, we agree with 



the results of reference [18| on this point.) 

This result should not be confined to strong coupling, but should persist 
for all coupling, i.e. the critical coupling in D = 1 should be at e 2 = 0. We 
conjecture that this sort of symmetry breaking for small e 2 is a manifestation 
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of the Peierls instability — the tendency of a one dimensional Fermi gas 
to from a gap at the Fermi surface. This happens with any infinitesimal 
interaction. 

In fact, this must also happen for the case where A^ is even. Then, there 
cannot be any spin order in 1 dimension. However, anomalies break the 
isoscalar chiral symmetry in the continuum theory and should also do so here. 
This means that there should be a dynamical generation of charge density 
wave which would be driven by the Peierls insability. It also implies that for 
Al = 2 for example, the ground state in the strong coupling limit would not 
be a Heisenberg antiferromagnet, but would be alternating empty site and 
site with two electrons in a spin singlet state. This state, even though it has 
large coulomb energy, avoids the infrared divergences of gapless Fermions. 

In higher dimensions, D > 2, it would be interesting to explore the pos- 
sibility of phase transitions between different symmetry breaking patterns 
for the SU(Al) flavor symmetry as one varies A"l- There is already some 
work on this subject in the condensed matter physics literature on SU(Al) 
antiferromagnets ]59]| . They analyze the SU(A^l) antiferromagnet which is 
similar to the strong coupling limit of an XJ(Nc) gauge theory (see [15J, |l6j for 
details) which is in the representation corresponding to a rectangular Young 
Tableau with Al rows and Ac columns. They work in the large Ac limit and 
show that there is a phase transition from the spin ordered Neel phase to a 
disordered phase when Al ~ Ac- In our case Ac = 1 so their analysis is not 
accurate. Nevertheless, we expect that there should be a phase transition to 
a disordered phase as Ac is increased. For odd A^ the chiral symmetry is 
always broken and the question we are asking is whether the flavor symme- 
try is also broken. For even Al possible phase transition is relevant to both 
chiral and flavor symmetry breaking. 

First of all, when A^c = 2, we have the j = 1/2 Heisenberg antiferromag- 
net which is known to have an ordered ground state in D > 3 and is also 
very likely to have an ordered ground state in D = 2. Furthermore, when 



Al — > oo, the ground state is known to be disordered in D = 2 [34] and is 
likely also the case in D = 3. In between there Al = 2 and Al = oo there 
should be a phase transition. It is interesting to speculate that the order- 
disorder phase transition which occurs as one increases Al in the SU(A" L ) 
antiferromagnet is the same one that appears in the study of chiral sym- 
metry breaking in 2+1-dimensional QED in the continuum [l^, |ll| where 
they find that chiral symmetry is broken only if the number of flavors is less 
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than a certain critical value. As we noted in the introduction, their work is 
effectively in the strong coupling (large e 2 ) limit. 

Our results also indicate that, besides the critical iV L , for a fixed N L 
which is small enough, there should also be a critical coupling constant e 2 
and, in fact, a critical line in the N^-e 2 plane where there is a second order 
phase transition between a spin ordered chiral symmetry breaking phase and 
a disordered (and possibly chirally symmetric phase). We speculate that in 
3+1-dimensions a similar situation could occur. 

A Notation 

In this paper we consider the lattice regularization most suited to the Hamil- 
tonian formalism where time is continuous and space is a square lattice with 
lattice spacing one. We use a finite spatial lattice Vr with lattice sites 

x = (xi, . . . , xd), —R < Xi < R (Al) 

were R is a positive integer and \Vr\ — (2R) D is the total number of lattice 
sites. In the thermodynamic limit, R — > oo. 

The lattice sites are connected by unit vectors 

1 = (1,0,...) 2 = (0,1,...) ... (A2) 

and the oriented link between the lattice site x and x + i is denoted [x, i]. The 
link oriented in the opposite direction is denoted —[x, i}. On the finite lattice 
Vr we identify lattice sites x and x + 2Ri. Then links obey the identity 

[x, — i] — — [x — i, i] (A3) 

The boundary of the [x, i] are the two points 

5[x,i] = (x + i) - x (A4) 

Also, a plaquette with corners x, x + i, x + i+j, x+j and with sides + 
[x + i + j, —i], [x + j, —j] is denoted as [x, and has the boundary 

5[x, i, j] = [x, i\ - [x, j] + [x + i, j] -[x + j, i\ ( A5) 
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It also obeys the identities 

= ~[x,i,j] [x,-i,j] = -[x-i,i,j] (A6) 

It is also possible to introduce higher dimensional structures, elementary 
cubes, etc. 

We shall also introduce lattice derivative operators, the forward difference 
operator 

Vif(x)=f(x + l)-f(x) (A7) 
and the backward difference operator 

VJ(x) = f(x)-f(x-l) (A8) 
The lattice Laplacian is 

V • V/(s) = E (f(x + i)- 2 f{x) + f{x - <)) (A9) 
i=i 

The functions which we shall consider are functions from either the lattice 
sites, links or plaquettes to the real numbers. The Fourier transform of a 
function on lattice sites is given by 

f(k) = J2 e lk - x f{x) (A10) 
xev R 

with / a function on the reciprocal lattice (momentum space) Vr with lattice 
sites 

k = (k±, . . . , k D ), ki = — — x integers, — n < fcj < n (All) 

2R 

The inverse Fourier transform is 



/(*)= I^Ee~*V(*) (A12) 



(note the the number of sites of the lattice Vr and its reciprocal lattice Vr 
are equal). The periodic delta functions on Vr and Vr are given by 

*(*)= l^lEe - ** (A13) 

Vr 



R\ 
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and 



5(k) = E (AM) 
key* 

and the Parseval relation is 

E r{x)g{x) = ±- E /*(-%(*) (A15) 

In momentum space, the lattice Laplacian is diagonal, 

(V^/)(*0 = -4Esin 2 (^/2)/(A;) (A16) 

i=l 

From this it follows that its inverse =^ is unambigously defined only on 

functions / with /(0) = J2 x ev R f( x ) = 0- We can extend its definition to all 
functions / by setting 

/(0) = (A17) 



-V - V 

>ace, the integra 

the Laplacian 



In position space, the integral kernel of _^ ^ is just the Green function for 



:x\ 



■\y)= E d Hx - v) T^D . m ^ ( A18 ) 



-V-V'^ fee ^ {0} 4Ef =1 sin 2 (^/2) 



For the infinite lattice V^, the momentum space is no longer a lattice but 
the Brillouin zone 

Voo = = {k = (ki, . . . , ko)\ki real , — n < ki < n} (A19) 

The inverse of the Fourier transform (|A10|) is then 

f(x) = I S^e-^f(k) (A20) 

JU B \2>'K) 

corresponding to the delta function on the infinite lattice, 

d D k 



s ^ = L^ e>kx (A21) 
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The Green function for the Laplacian on Voo is given by 

wz^> - L-^^m^m (A22) 

and is well defined in dimensions D > 2. In two dimensions it is defined 
by one additional subtraction which removes the logarithmic divergence in 
the integration. In one dimension it must be defined by solving the Laplace 
equation with a source explicitly. 

In the main text we appreviate J2xeV R as an d similarly in momentum 
space. 



B Lattice Fermions 

Throughout this paper we shall use the staggered Fermion formalism which 
was originally developed by Kogut and Susskind. This formalism is well 
known and the details can be found in the papers of Susskind and collab- 



orators ||0|, [L7|, [1^, [L9| and Kluberg-Stern et.al. |4l]]. Here we shall review 
the basic features and make some observations which are necessary for our 
present discussion. Some of these observations have already been made in 



We shall use the staggered Fermion formalism since we believe that it 
gives the closest possible analog to the lattice Fermions encountered in con- 
densed matter physics. As a regularization of Fermions in relativistic quan- 
tum field theory, this formalism has the disadvantage that chiral symmetries 
are discrete, rather than continuous. The method should be regarded as 
adding some formally irrelevant operators to the Hamiltonian. These opera- 
tors make the Hamiltonian local but break the continuous chiral symmetry 
down to a discrete subgroup. (Actually, there is a non-local chiral symmetry. 
However, being non-local it is not a useful symmetry in that, for example, 
it does not imply the existence of Goldstone Bosons in the phase where it is 
broken.) 

Thus, we can really only address questions about discrete chiral symmetry 
breaking. This should be enough to tell us whether mass generation, and in 
fact what sort of mass generation, is possible. 
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B.l Review of Staggered Fermions 



The purpose of the staggered Fermion method is to minimize Fermion dou- 
bling which always accompanies lattice Fermions. Generally, staggered Fermions 
are obtained by the spin-diagonalization method. To implement this method, 
we begin with the naively latticized Dirac Hamiltonian, 

H F = - (ip\x)ia j V 'jip(x) - (V^X^'o^x)) (Bl) 

where a J are the 2 tO * 1 )/ 2 ! -dimensional Dirac a-matrices. (Here [(D + l)/2] 

is the integer part of (D + l)/2.) They are Hermitean, = a- 7 and obey 
the Clifford algebra 

{a\a j } = 25 i3 (B2) 

They are therefore unitary matrices, a 1 ' a? = 1. Using the properties of the 
difference operator, flBTf) can be presented in the form 

R f = ~\ £ (f^ + 5V>0*0 - ^{x)a^{x + j)) (B3) 

Since the Dirac matrices are unitary, the naive lattice Fermion Hamiltonian 
in ( P3|) resembles a condensed matter Fermion hopping problem with a back- 
ground U(2^ D+1 ^ 2 ^) gauge field given by the a-matrices. In any plaquette of 
the lattice, [x,i,j], this background field has curvature 

c/aW't = -1 (B4) 

The curvature resides in the U(l) subgroup of C/(2^ £)+1 ^ 2 I) and has exactly 
half of a U(l) flux quantum per plaquette. This is true in any dimensions. 
We observe that either 1/2 or zero flux quanta are the only ones allowed 
by translation invariance and parity and time reversal symmetries of the 
Hamiltonian. 

Since the curvature if the a- matrices is U(l)-valued, we should be able to 
do a gauge transform which presents the matrices themselves as U(l) valued 
gauge fields (i.e. diagonal). A specific example of such a gauge transform 
due to Kluberg-Stern et. al. [41| is 

V>(x) -> (a 1 )" 1 (a 2 )" 2 . . . (a D ) XD ij>(x) (B5) 
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Then 



^(x + j)a j i[)(x) -> (-1)^=1 x ^(x + j)i[)(x) (B6) 
The resulting Hamiltonian is 



R f = -^E(- 1 ) ECllXfc (tf(x+jMx)-tf(xMx + jj) (B7) 



This describes 2^" D+1 ^' 2 ' identical copies of Fermions with the same Hamil- 
tonian which must all give Fermions with the same spectrum as the original 
Hamiltonian in ( |B1| ). Staggered Fermions are obtained by choosing one of 
these copies. This reduces the Fermion doubling by a factor of the dimension 
of the Dirac matrices, 2^ D+1 >' 2 K 

In the staggered Fermion method, we treat the components of the original 
lattice Dirac Hamiltonian as flavors, rather than components of the relativis- 
ts spinor necessary for Lorentz invariance. The spinor components now 
reside on adjacent lattice sites. In this method, the continuous chiral sym- 
metry of the massless Hamiltonian, under the transformation tp — > e* 7 9 ip 1 is 
lost. There is a discrete chiral symmetry, corresponding to translations by 
one lattice site in any direction. Explicitly, 

ip(x) -> (-l^^+^VO + j) (B8) 



is a symmetry of the Hamiltonian (|B7|) and corresponds to a discrete chiral 
transformation. 

Mass operators correspond to staggered charge densities. The operator 
E = 53(-l)Sfc=i**V t (^)V'(aj) (B9) 



changes sign under the chiral transformations (|B8|) and corresponds to a 
certain Dirac mass. 

With staggered Fermions there is still a certain amount of Fermion dou- 
bling. The doubling can be counted by noting that the staggered Fermion 
Hamiltonian (|B7|) is invariant under translations by two lattice sites. There- 
fore, a unit cell is a unit hypercube of the lattice, containing 2 D sites and 
staggered Fermions correspond to a 2 D component spinor. The dimension 
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of the Dirac matrices is d^ D+1 '' 2 *. Therefore the number of Dirac spinors 
we obtain is 2 D /2^ D+l ^ 2 \ For lower dimensions the minimum number of 
continuum flavors can be tabulated as 

d dim. of Dirac matrices No. of flavors 

1 + 1 2 1 

2 + 1 2 2 

3 + 1 4 2 

Only in 1+1- dimensions do we get a single species of Dirac Fermion. 

B.2 Explicit Example in 3+1 Dimensions 

For simplicity in notation, the formulas here and in the following subsection 
are given for an infinite spatial lattice V^. 

To see how to take the continuum limit explicitly, consider the case of 
d = 3 + 1. There, we divide the lattice into eight sublattices and label the 
spinor components as 

'(/'(even, even, even) = ip 1 ^(odd, even, odd) = ip 7 (BIO) 

ip (even, odd, even) = ip 6 ^(even, even, odd) = ip§ (Bll) 

ijj (odd, odd, even) = ip A ^(odd, even, odd) = ip 3 (B12) 

ijj (even, odd, odd) = ip 2 (odd, odd, odd) = ip$ (B13) 
In terms of these spinors, the Hamiltonian (|B7[) can be written as the matrix 



operator 

H= f d 3 kt/j J( (k)A i sin kii>(k) (B14) 

J fig 

where 

Q B = {ki : -tt/2 < ki < vr/2} (B15) 
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is the Brillouin zone of the (even,even,even) sublattice, 

' a' 
a 1 



^=L°* n) (B16) 



and 



are a particular representation of the Dirac matrices. 
In this representation the mass operator is 

d 3 k 



B#=>V(#)= —-^{k)B^{k) (B18) 

Jn B {Zn)° 



where 

-1 



fi=(i ° ) (B19) 



The Fermion spectrum is 
uj{k) : 



. sin 2 1* + m 2 (B20) 



and only the region ki ~ is relevant to the continuum limit. We have 
normalized ip{k) so that 



{^(x),^(y)}=6(x-y) , {^(k), ^(l)} = 8(k - I) (B21) 



If we define 



■2 







0=( O ->) (B22) 



M=U!7^ ! + ^) (B23) 



and the unitary matrix 

M — 
and 

i> = M^)' (B24) 
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with 



tf = U> a ,iM (B25) 

the Hamiltonian is 

"<=L*wM''* 1 *- pm J (*) < B26 > 

In the low momentum limit, sinfcj ~ hi, with Fermion density 1/2 per site so 
that the Fermi level is at the intersection point of the positive and negative 
energy bands, we obtain 2 continuum Dirac Fermions. 

This describes two flavors of 4-component Dirac Fermions and the Dirac 
masses for each component given by the staggered charge density have op- 
posite signs. Thus the charge density breaks the discrete chiral symmetry. It 
also breaks a flavor symmetry which, in the absence of mass, mixes the two 
continuum Fermions. 



B.3 General Continuum Limit 

We shall now consider the continuum limit in a general number of dimensions. 
A formalism much like (but not exactly the same as) the present one can be 
found in [f4l . 



We shall begin with the Hamiltonian 



We consider an elementary hypercube of the lattice with sides of length 1 
and 2 d sites generated by taking a site all of whose coordinates are even and 
adding to it the vectors 

a = {a u ...,a D ) on = OR 1 (B27) 

We also decompose the lattice into 2 D sublattices generated by taking a site 
of the elementary hypercube, (even,even,...)+<3 for some a and translating it 
by all even multiples of lattice unit vectors, i. We label the Fermions which 
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reside on the sublattice of each of the corners of the elementary hypercube 
as ip ai ,...,a D {x)- In momentum space the Hamiltonian is 



H f = f d D kJ2 ^ ai _ aD (k) I';,,..,,,, sinfc, ipp x ...p D {k) (B28) 



where Qb = {ki : — 7r/2 < ki < vr/2} is the Brillouin zone of the (even, even,. . .) 
sublattice, the momentum space Fermions have the anticommutator 

{^...a D {k)^l^ D {k')} = (n<W,) S(k - k') (B29) 

and the Dirac tensors are 

r L 1 ...a D /3 1 .../3 = <Wi • • • ^.-ift-j^ft^+ift+i • • • S <*dPd (-l)^*=l °* (B30) 

They obey the Clifford algebra 

FT' + I'M'' = 25 ij (B31) 



The spectrum of the Dirac operator is u(k) = iyZ^Li sin 2 ki. Note that, to 
set the Fermi level of the Fermions at the degeneracy point where the two 
branches of the spectrum meet, it is necessary that the Fermion states are 
exactly half-filled. This is also required for charge-conjugation invariance, or 
particle-hole symmetry of the vacuum state. 

The staggered charge density operator ( jB9|) is equivalent to a mass oper- 
ator where 

£= / d J, ArVi 1 „^(A ; )r^„. ai , A „^Vft..^(fc) (B32) 

where 

r°,..a Dft ... to = <W, • • •S anPn {-ljZZ*'"> (B33) 

Here, T satisfies the algebra 

r o r i + T i T o = q ( B34 ) 

and 

r°r° = 1 (B35) 
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Thus the spectrum of the operator Hf + mE is 



u(k) = ± 



\ 



D 



sin h 



(B36) 



i=l 



which is the spectrum of a relativistic Fermion in the limit k ~ 0. 

Here, we count the number of flavors of Fermions obtained in the con- 
tinuum limit by noting that (|B28|) describes a 2 jD -component Fermion. In 
D dimensions the Dirac matrices are [(D + l)/2] dimensional, therefore the 
continuum limit of (|B28|) describes 2 D /2^ D+1 ^ 2 ^ species of Dirac Fermions. 
These are tabulated up to dimension 4 in (OTl). 



C Fermion Representation of SU(iV) Quan- 
tum Antiferromagnet 

The Hamiltonian for an U(iV) quantum antiferromagnet is 

g 2 

H A FM = — J ab{x)Jba{y) (CI) 

<x,y> 

where Jab(x), a, b = 1, • • • , iVx, obey current algebra relations associated with 
the Lie algebra of XJ(N), 

[Jab(x), Jcd{y)] = (J a d(x)5 bc - J cb (x)5 ad )S(x - y) (C2) 

and where < x, y > denotes the link connecting sites x and y on a bipartite 
lattice. For simplicity, we shall take the lattice to be cubic. Here, we have 
used a particular basis for the SU(iV) algebra which can be conveniently 
represented by Fermion bilinear operators, 

J ab {x) = r\x)^\x)-5 ab /2 (C3) 

The representation of the algebra on each site x is fixed by specifying the 
Fermion number of the states, 

P (x)=J2Jaa(x) (C4) 
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For example, the Fermion vacuum state |0 > such that 

■^ a (x)\0 >= , Va,x (C5) 
is the singlet state, the states with m < N Fermions per site, 
[]^ ait (a;)^ a2t (a;)...^ amt (^)|0 > 

X 

corresponds to p{x) = m — N/2 for all x and the irreducible representation 
with the Young Tableau 



Figure 3:When there are m Fermions per site, the representation of SU (Nl) 
has Young Tableau with one column of m boxes. 

For each site x, p(x) is the generator of the U(l) subgroup of U(iV). Using 
a basis T l = (T*)*, i = 1, . . . , N 2 — 1, of the Lie algebra of SU(iV) in the 
fundamental representation normalized so that ti{T l T^) = T^ b T^ a = 5^/2, 
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and using 



T ab T cd = ^adSbc - J^$ab$cd (C6) 

it is convenient to introduce 

j\x) = r\x)r ah r{x) (C7) 

obeying current algebra of the Lie algebra of SU(iV), and to write the Hamil- 
tonian (CT) as 



Hafm = t; E P{x)p{y) + Hsv(N) (C8) 
iV <x,y> 

with 

H S v(N)=9 2 £ JK^JKv) (C9) 

<x,y> 

is the Hamiltonian of an SU(iV) antiferromagnet. From this it is obvious 
that by fixing the p(x), Hafm is reduced to an SU(A^) antiferromagnet. 

For example, the familiar j = 1/2 SU(2) Heisenberg antiferromagnet is 
obtained from ( PH) , = 2, by using the identity 



corresponding to ( C6[) for N = 2. 



Generally, when iV is even we will consider the representations where 
m = N/2, so that the Fermion occupation of each site is N/2 and p(x) = 0. 
When iV is odd we divide the lattice into two sublattices such that the 
nearest neighbors of all sites of one sublattice are in the other sublattice 
(when this is possible the lattice is said to be bipartite). When iV is odd, 
the representation of SU(iV) has (iV+ l)/2 Fermions, i.e. p(x) = 1/2, on the 
sites of one sublattice and (iV — l)/2 Fermions, i.e. p(x) = —1/2 on the sites 
of the other sublattice. 
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